Superfluid transition in a correlated defect network 
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Motivated by recent experiments on possible supersolid behavior of 4 He solids at low temperature, 
we consider a model of superfluidity in a defected solid containing a system spaning network of 
correlated linear dislocations, or planar grain boundaries. Using arguments based on the Harris 
criterion, as well as numerical simulations, we find that such correlated quenched disorder shifts the 
familiar superfluid lambda transition to a new disordered universality class in which the correlation 
length exponent v > 1. This results in the temperature-derivates for the superfluid density, dp s /dT, 
and for the heat capacity, dc/dT, remaining finite at the transition T c , and thus a less singular 
transition, profoundly different from the usual lambda transition. 



PACS numbers: 64.60.A-, 67.80.bd, 67.25.dj 

Considerable excitement has been generated by the ob- 
servation by Kim and Chan (TJ [2] , that the period of os- 
cillation of a torsional oscillator (TO) containing solid 
4 He decreases as the system is cooled below roughly 250 
mK. Two different scenarios have been suggested to con- 
tribute to this phenomenon, (i) Initially, it was suggested 
that a transition to a supersolid state resulted in a reduc- 
tion of the mass viscously coupled to the oscillator walls 
and hence a nonclassical rotational inertia (NCRI). Sub- 
sequent experiments showed this NCRI was closely tied 
to the presence of defects in the 4 He crystal: the slower 
the crystal was annealed, the smaller was the resulting 
NCRI signal [3]. Numerous theoretical calculations and 
simulations followed, supporting the idea that a super- 
solid should not exist for a pure 4 He crystalline state 
0HB], however superfluidity could exist in the cores of 
crystalline defects [8] . A system spanning network of 
pinned defects could then lead to an effective supersolid. 
(ii) Alternatively, it has been suggested that the appar- 
ent NCRI is a classical effect due to the viscoelasticity of 
the defected 4 He solid [HI [TO] • At high temperatures de- 
fects are mobile and a fraction of the 4 He lags the driving 
oscillation. But upon cooling, a glassy state is entered in 
which the defect relaxation time r grows dramatically. 
When the relaxation rate 1/r decreases below the TO 
driving frequency w, the defects lock into phase with the 
driving oscillation, leading to a shift in the TO period. 

While recent experimental works [TTJ [TO] seem to favor 
(ii) as the dominant mechanism for the observed NCRI, 
this scenario still assumes a highly defected crystal; one 
might therefore still expect to see some contribution from 
superfluidity in the defect network, as in scenario (i). 
While direct probes of a steady state superfluid flow have 
yielded conflicting results [TO], other experiments have 
failed to see the characteristic features expected for the 
bulk lambda transition of three dimensional (d = 3) su- 
perfluid 4 He [H]. In particular, (1) the TO period shift 
varies rather smoothly with temperature [U El [T51 [16] , 
with no sign of the sharp singularity expected if some part 



of this shift was due to a superfluid density p s ~ \T—T C \ V , 
with v w 0.67; and (2) equilibrium measurements of heat 
capacity [IT] [TO] show a smooth bump rather than the 
sharp lambda-shaped cusp characteristic of a bulk super- 
fluid transition with c ~ \T - T c \~ a , a = 2 - dv « -0.01 
[19] , One possible reason for the absence of such sharp 
features is if the dislocation density is simply too small 
to give an observable superfluid signal [7 . In this work, 
however, we suggest another possible contributing fac- 
tor: that the form of the disorder may shift the super- 
fluid transition to a new universality class, with much 
less singular behavior. 

The Harris criterion [30] argues that uncorrelated 
point-like disorder is irrelevant whenever 2 < dv, with 
d the spatial dimension of the system, and v the corre- 
lation length critical exponent. For bulk superfluidity in 
d = 3, the pure system has f purc ~ 0.67 and the Harris 
criterion is satisfied. It has therefore been argued [2"Tl - f2"3] 
that the superfluid onset in a random dislocation network 
should show the same lambda singularity as in bulk su- 
perfluid helium. However, since the dislocation cores are 
continuous one dimensional objects, the assumption of 
uncorrelated point-like disorder might not be appropri- 
ate. In this work we consider the case when disorder is 
correlated along a network of intersecting lines or planes. 
In all cases we argue that the effects of such correlated 
disorder is to lessen the sharpness of the critical singu- 
larity, removing the divergent temperature-derivatives at 
T c that are characteristic of the usual bulk lambda tran- 
sition. 

To model superfluidity in the dislocation cores of a de- 
fected 4 He crystal, we start with a lattice 3D XY model. 
Neglecting amplitude fluctuations of the condensate or- 
der parameter 4 r (r) = |^|e !e W, the fluctuations in the 
phase angle 9 that drive the superfluid transition can be 
modeled by the Hamiltonian, 

W = -5^J |i (r i )cos[fl(r i + A)-«(r0] , (1) 
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where i labels the discrete sites of a simple cubic lattice 
of length L with periodic boundary conditions in all di- 
rections, and J M (i"i) is the coupling on the bond connect- 
ing site Vi to Yi + /2, with /i = x,y, z the lattice axis direc- 
tions. We will choose spatially inhomogeneous random 
couplings J^(rj) according to several different schemes. 
Although in real solid 4 He dislocation cores can thermally 
fluctuate at higher temperatures, we will assume that at 
sufficiently low T the dislocations are pinned by 3 He im- 
purities, and so we take the couplings J„ (rj) to be fixed 
quenched variables. 

According to Harris [20] , the effect of disorder is con- 
trolled by the mean square fluctuation of the local cou- 
pling A J, disorder averaged over a sub-volume V = R d 
of the system, 



[(AJ)|] = 



R 



, (2) 



where [. . . ] denotes the average over different realizations 
of quenched disorder. For a continuous transition with 
diverging correlation length £ ~ \T — T c \~ v , the Harris 
criterion argues that disorder is irrelevant at T c when the 
coupling fluctuation averaged over a correlation volume 

£ rf satisfies y^[(AJ)|] < \T - T c \ as T -> T c . For future 

use, we note that [(A J)^] can be rewritten in terms of 
the coupling correlation averaged over the sub- volume V, 

l(Ajf R ] = C(R) ^R- d J2J2 [SUrdSJ^m , (3) 

At, At' Ti£V 

where 5J M (rj) = J M (rj) — [J^J. For uncorrelated point dis- 
order, with [<5J M (ri)5J^'(0)] oc <5 MiM '(5 rii0 , Eq. ([3} yields 
C(£) - ~ |T - T c \ dlJ . The Harris criterion for 

the irrelevance of such point disorder then becomes 
\T - T c \ dy ' 2 <\T- T c \ or the familiar 2 < dv. 

We now wish to consider models of correlated disor- 
der. One such case was considered many years ago in 
a seminal work by Weinrib and Halperin [21]. They 
noted that when spatial disorder is introduced as a set 
of straight randomly oriented lines, the resulting disor- 
der averaged coupling correlation decays algebraically as 
[6J{r)SJ(0)} ~ r~ a , with a = d-l. In this case, Eq. {3 
gives C(£) ~ £~ a and the Harris criterion argues that dis- 
order is irrelevant when 2 < av. Using a renormalization 
group expansion for weak Gaussian disorder, Weinrib and 
Halperin showed that this is indeed the case: when a < d, 
disorder is irrelevant whenever 2 < av puro . They further 
found that when this criterion fails, the correlation length 
exponent at the new disordered critical point satisfies 
^disor = 2/a. Applying this conclusion to a 3D supcrfluid, 
we see that such linear disorder with a = 2 is relevant, 
resulting in a superfluid density that vanishes linearly as 
p s ~ \T — T c | l ' diaor , with fdisor = 1) an d a heat capacity 
c ~ IT — T c \~ a with exponent a = 2 — elisor = — 1- In 



particular, dc/dT scales as \T— T c | _a_1 , with — a— 1 = 0, 
and thus remains finite at T c in contrast to its divergence 
in the pure model. Thus for such linearly correlated dis- 
order, we expect that the superfluid density vanishes less 
sharply than the pure model, while the heat capacity has 
a far less singular cusp. For planar disorder, as might be 
the case if superfluidity is carried on randomly oriented 
twin grain boundaries, a = d — 2 = lin 3D, and the 
superfluid singularities are even less sharp. 

We next consider a slightly different model of linearly 
correlated disorder. Motivated by the notion that dis- 
location lines in solid 4 He may preferentially align par- 
allel to crystalline axes, we consider a 3D XY model in 
which correlated disorder enters along straight lines ori- 
ented only along the three lattice directions. We choose 
quenched couplings as follows. Within the plane at r p = 
0, we randomly choose equal numbers of bonds J M from 
the bi-valued distribution J = I ±6. We then continue 
these couplings in correlated straight lines by requiring 
J^Ti + p,) = Jfj,(ri). We follow this procedure for all 
three directions /i = x,y, z. For the results presented be- 
low we use 5=1; superfluidity is strictly confined to the 
dislocation cores and the three dimensionality of the su- 
perflow results solely from the intersections of these cores 
to form an interconnected network. We have also consid- 
ered the case 5 = 0.95, as a model in which superfluid 
particles may tunnel through the bulk between disloca- 
tion cores. We find the critical behavior to be the same 
in both cases. The coupling correlation for this disorder 
is [SJ^SJ^O)] cx [5(x)5(y)+S(y)5(z) + S(z)S(x)]d litfl >. 
However the volume averaged correlation is C(R) ~ R~ a 
with a = d — 1, just as for randomly oriented lines. We 
thus expect this disorder to be relevant; applying the 
Harris criterion at the new disordered critical point then 
requires fdisor > 2/a = 1. 

To confirm this behavior, we carry out extensive nu- 
merical simulations. The superfluid density is propor- 
tional to the XY helicity modulus, which (in units where 
m/fi = 1) is given by [25] 



Ps = L 
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J2 J z (r i )cos[e(r l + z) - 0(n 
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The XY magnetization is M = |£ l e lf,(r,) | 
Binder cumulant U is given by, 

'{MY 

We also consider the heat capacity c, 

c = L -3 T -2 [ (h2) _ (h) 2] 

and the XY spin susceptibility x, 

X = L^T- 1 [(M 2 ) - (A/) 2 ] 



(4) 



and its 



(5) 

(6) 
(7) 
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To analyze the critically of the transition at T c we ap- 
ply finite-size-scaling (FSS) methods [26]. We expect p s , 
U, and x as a function of T and system length L, to obey 
the usual FSS relations, 



Ps (T, L) = L^nttL 1 '") , U(T, L) = UitL 1 ^) 
X {T,L)=Lil»X{tL 1 l») , 



(8) 



where IZ(-), U(-) and X(-) are scaling functions, and t = 
(T-T c )/T c . 

To carry out our Monte Carlo (MC) simulations to high 
accuracy we used an effective Wolff collective update al- 
gorithm, which minimizes the effects of critical slowing 
down at the transition [57] . As the Wolff algorithm is not 
so efficient at lower temperatures, we add one ordinary 
Metropolis sweep through the system for each 100 Wolff 
updates. Data for the averages were accumulated dur- 
ing 4 - 16 x 10 3 Monte Carlo steps (MCS), after equally 
many initial MCS were discarded to reach equilibrium. 
One MCS is defined as L 3 single site update attempts. 
We tested for equilibration by increasing the number of 
discarded initial MCS until stable results were obtained. 
For our FSS analysis, quenched disorder averages were 
computed over 10 3 — 10 4 independent realizations of the 
random couplings J^r^). 

The main qualitative results of this work are shown in 
Fig. [I] where we plot the superfluid density p s and heat 
capacity c vs T/T c of a system of fixed length L = 40, 
for both the pure 3D XY model (T c = 2.203) and for 
our model of Eq. (JlJ with linearly correlated disorder 
(T c = 2.501, as determined below). Results for the disor- 
dered case are averaged over several hundred realizations 
of the quenched random couplings. We see that the pres- 
ence of the disorder dramatically softens the singulari- 
ties of the pure system, removing the sharp singularities 
of the lambda transition. Derivatives with respect to 
temperature that diverge as T — » T c in the pure model, 
appear to become finite in the disordered model. 

To verify this conclusion, and to demonstrate that the 
model with linearly correlated defects belongs to a new 
disordered universality class, we perform a FSS analysis 
to determine critical exponents [35]. In Figs.[2ja) and (c) 
we plot respectively Lp s and U vs T for different system 
sizes L = 10 — 80. Equations (JsJ) predict that the curves 
for different L should all intersect at the common point 
t — 0, i.e. when T = T c . We see that as L increases, the 
curves do indeed seem to be approaching a common in- 
tersection point, yielding T c « 2.5. The deviations from 
a perfect common intersection are due to corrections to 
scaling, which can be noticeable when L is insufficiently 
large [19] . 

For a more accurate determination of the critical T c , 
as well as the correlation length exponent v, we expand 
the scaling functions of Eqs. ^ as third order polynomi- 
als, and fit our data to these scaling forms with T c , v and 
the polynomial coefficients as free fitting parameters. We 
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FIG. 1. (color online) Superfluid density p a (solid symbols) 
and heat capacity c (open symbols), vs T/T c , for several 
model systems on a lattice of length L — 40. Black circles: 
pure 3D XY model; red squares: straight line disorder; green 
diamonds: directed random line disorder; blue triangles: ran- 
dom planar disorder. For the disordered models, results are 
averaged over several hundred to more than a thousand ran- 
dom coupling configurations 
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FIG. 2. (color online) (a) Scaled superfluid density Lp s , (c) 
Binder cumulant U, and (e) susceptibility x vs T for system 
sizes L — 10 — 80. Scaling collapse of (b) Lp E , (d) U, and (f) 
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2.501, v = 1 and 7 = 1.84. 



tL}l u , for system sizes L — 30 — 80, using values 



restrict our fits to data satisfying \tL x l v \ < 0.5. Using 
system sizes L m - m to £ max = 80, we plot the resulting 
fitted values for T c and v vs 1/L min in Figs.[3^a) and (b) 
respectively. We see that as L m i n increases, the values of 
T c and v from Lp s and U approach each other. A pro- 
cedure including corrections to scaling yields consistent 
values of T c = 2.501 ± 0.001 and v = 1.00 ± 0.05 from 
fits to both Lp s and U. Errors represent one standard 
deviation statistical error as estimated using the method 
of synthetic data sets [55]. Using these common values 
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FIG. 3. (color online) Values of (a) T c and (b) v obtained 
by fitting data with \tL x ^\ < 0.5 to third order polynomial 
expansions of the scaling functions of Eqs. pjl, using data 
from system sizes L m j n to L max = 80. 



of T c and v, the resulting scaling collapses, plotting Lp s 
and U vs tL ■' v , are shown in Figs, gb) and (d), for sizes 
L = 30 — 80. We note that our result v s» 1 is identical to 
the prediction of Weinrib and Halperin for algebraically 
correlated couplings with a = d— 1 = 2, and implies that 
the new disordered fixed point just marginally satisfies 
the Harris criterion for stability. 

Next we consider the susceptibility \ of Eq. In 
Fig. [2je) we plot x vs T for system sizes L = 10 — 80. 
The exponent 7 in Eq. |8]) is obtained from a power law 
fit to the maximum value of x(T) vs L. We find 7 = 
1.84 ±0.1. Using this value of 7, and T c = 2.501, v = 1.0 
as determined above, we show in Fig. [21(f) the resulting 
scaling collapse, plotting xL^ 1 ^ vs tV^ 1 for system sizes 
L = 30 — 80. The collapse is excellent. From the two 
independent exponents v w 1-0,7 ~ 1-84 all the other 
thermodynamic critical exponents can be obtained from 
standard scaling laws [25]. Hyperscaling gives the heat 
capacity exponent a = 2 — dv k, —1.0; the Rushbrook 
equality, a + 2/3 + 7 = 2, gives for the order parameter 
exponent (3 ~ 0.58; and the correlation function exponent 
•q is obtained from 7 = (2 — rj)v giving r\ ss 0.16. 

So far, our linear disorder has taken the form of per- 
fectly straight lines. In a solid, however, dislocation cores 
may wander as they traverse the system. To consider the 
effect of such wandering, we generalize our model by let- 
ting each defect line be a directed random walk. For a 
walk directed along z, for example, each step in the z 
direction is allowed to include a random transverse fluc- 
tuation Ar_L, sampled equally from Arj_ = {0, ±x,±y}. 
For each walk the Arj_ are constrained to sum to zero, 
so that the line at z — L returns to its starting position 
at z — under the periodic boundary conditions. For 
each such walk j, a variable rij — 1 is placed on each 
bond of the walk. The couplings of the corresponding 
XY model are then set to J = J2j n j on each bond. Fi- 
nally, all couplings J are rescaled by a constant factor 
so that the disorder average [J M ] = 1. We thus con- 
struct such configurations on a lattice of size L = 200 and 



numerically compute the integrated coupling correlation 
C{R) of Eq. (J3j>, averaging over more than 1000 differ- 
ent realizations of the disorder. Our results are shown 
in Fig. [3j along with C(R) as computed for our original 
straight line model. We find a clear algebraic decay in 
both cases, with C(R) ~ R~ a , a = 2. Thus the directed 
random line model is expected to be in the same univer- 
sality class as our straight line model. In Fig. [T] we show 
results for p s and heat capacity c for this directed random 
line model, as obtained from simulations on a lattice of 
size L = 40, averaging over several hundreds of disorder 
configurations. The effect of superfluidity seems slightly 
suppressed, as compared to the straight line case, but the 
shapes of the singularities appear to be the same. 
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FIG. 4. (color online) Integrated coupling correlation C(R) 
of Eq. j3j, vs R, for a system of length L = 200. Results 
are normalized by C(l), and shown for the case of straight 
randomly positioned lines, directed random lines, and random 
flat planes. 



We have also considered a random plane model. We 
select a sequence of flat planes as follows. For planes ori- 
ented with normal along z, each xy plane at height z is 
randomly selected or not with probability 1/2. Such ran- 
dom plane sequences are selected in all three directions x, 
y, z simultaneously. For each such plane j we set nj = 1 
for all bonds in the plane; the XY bond couplings are 
J = Y^j n j, an d then rescaled so that [J] = 1. In Fig. [4] 
we plot the resulting numerically computed correlation 
C{R). We find an algebraic decay with C(R) ~ R~ a , 
0=1. In Fig. [T] we show results for p s and heat capacity 
c for this random plane model, as obtained from simu- 
lations on a lattice of size L = 40, averaging over more 
than one thousand disorder configurations. As expected, 
the singularity at T c is now even smoother than for the 
random line models. 

To conclude, we have considered a variety of models for 
a superfluid with long range correlated quenched disor- 
der. In all cases we have argued that the disorder changes 
the universality class of the superfluid transition to one 
with v > 1. Our results show that superfluidity in such 
correlated disorder networks, such as have been proposed 
in some models for solid 4 He, would not display the famil- 
iar sharp features of the lambda transition of an ordinary 
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bulk supcrfluid. 

We are grateful to Egor Babaev, Alexander Bal- 
atsky, and John Reppy for discussions. This work 
was supported by the Swedish Research Council, the 
Goran Gustafsson foundation, NSF grant DMR-1205800 
and the Swedish National Infrastructure for Computing 
(SNIC). 



[1] E. Kim and M. H. W. Chan, Science 305, 1941 (2004). 
[2] E. Kim and M. H. W. Chan, Nature 427, 225 (2004). 
[3] A. S. C. Rittner and J. D. Reppy, Phys. Rev. Lett. 98, 

175302 (2007). 
[4] N. V. Prokof'ev, Advances in Physics 56, 381 (2007). 
[5] B. K. Clark and D. M. Ceperley, Phys. Rev. Lett. 96, 

105302 (2006). 

[6] M. Boninsegni, N. V. Prokof'ev, and B. V. Svistunov, 

Phys. Rev. Lett. 96, 105301 (2006). 
[7] M. Boninsegni, A. B. Kuklov, L. Pollet, N. V. Prokof'ev, 

B. V. Svistunov, and M. Troyer, Phys. Rev. Lett. 99, 

035301 (2007). 

[8] S. G. Soyler, A. B. Kuklov, L. Pollet, N. V. Prokof'ev, 
and B. V. Svistunov, Phys. Rev. Lett. 103, 175301 (2009). 
[9] J. Day and J. Beamish, Nature 450, 853 (2007). 
[10] A. V. Balatsky, M. J. Graf, Z. Nussinov, and J.-J. Su, 

arXiv:1209.0803. 
[11] J. D. Reppy, Phys. Rev. Lett. 104, 255301 (2010). 
[12] A. D. Fefferman, X. Rojas, A. Haziot, S. E. Balibar, 
J. T. West, and M. H. W. Chan, Phys. Rev. B85, 094103 
(2012). 



[13] M. W. Ray and R. B. Hallock, Phys. Rev. Lett. 100, 
235301 (2008). 

[14] M. Barmatz, I. Hahn, J. A. Lipa, and R. V. Duncan, 

Rev. Mod. Phys.79,1 (2007). 
[15] J. T. West, O. Syshchenko, J. Beamish, and M. H. W. 

Chan, Nature Physics 5, 598 (2009). 
[16] A. C. Clark, J. T. West, and M. H. W. Chan, Phys. Rev. 

Lett. 99, 135302 (2007). 
[17] X. Lin, M. H. W. Chan, and A. C. Clark, Nature 449, 

1025 (2007). 

[18] X. Lin, A. C. Clark, Z. G. Cheng, and M. H. W. Chan, 

Phys. Rev. Lett. 102, 125302 (2009). 
[19] M. Campostrini, M. Hasenbusch, A. Pelissetto, and E. 

Vicari, Phys. Rev. B74, 144506 (2006). 
[20] A. Harris, J. Phys. C: Solid State Phys. 7, 1671 (1974). 
[21] J. Toner, Phys. Rev. Lett. 100, 035302 (2008). 
[22] D. Goswami, K. Dasbiswas, C.-D. Yoo, and A. T. Dorsey, 

Phys. Rev. B 84, 054523 (2011). 
[23] M. Boninsegni and N. V. Prokof'ev, Rev. Mod. Phys. 84, 

759 (2012). 

[24] A. Weinrib and B. I. Halperin, Phys. Rev. B27, 413 
(1983). 

[25] Y. H. Li and S. Teitel, Phys. Rev. B 41, 11388 (1990). 

[26] V. Privman, in Finite Scaling and Numerical Simulations 
of Statistical Systems, edited by V. Privman (World Sci- 
entific, Singapore, 1990). 

[27] U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

[28] W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. 
Flannery, Numerical Recipes 3rd ed. (Cambridge Univer- 
sity Press, 2007). 

[29] S.-K. Ma, Modern Theory of Critical Phenomena (Fron- 
tiers in Physics, 1976). 



